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MATRIX  DATA  ANALYSIS 
COLOR/B&W  CODING  IS  NOT  ALWAYS  ENOUGH 


Background 


<5 


The  Defense  Mapping  Agency  produces  digital  data  bases  that  describe  the 
physical  appearance  of  the  surface  of  the  earth.  These  data  bases  include,  but 
are  not  limited  to,  terrain  elevation,  culture  including  landscape  character¬ 
istics,  and  vertical  features.  This  data  is  collected  from  digitized  source 
maps,  from  optically  or  digitally  correlated  stereo-pairs  of  photographic 
imagery,  and  from  digital  multi-spectral  sensor  data.  A  dramatic  impact  has  been 
made  in  the  ability  to  analyze  these  digital  data  bases  by  applying  state-of-the- 
art  digital  image  technology  processing  and  display  concepts.  These  include  a 
variety  of  color  and/or  black'  and  white  displays  of  not  only  intensity/color 
coded  matrix  data,  but  also  image  processed  data  using  specialized  convolution 
filters,  texture  discrimination,  and  special  color  representation  techniques. 
In  addition,  computer  generated  imagery  from  these  data  bases  serves  as  a  final 
analysis  tool. 

A 

Techniques  andVResults 


With  results  similar  to  those  from  manually  compiled  classical  mapping  and 
charting  displays,  digital  terrain  elevation  data  may  be  used  to  generate  a 
standard  contour  plot  (figure  1.)  and  the  corresponding  tint  plate  (figure  2.)  in 
which  the  areas  between  contour  lines  are  either  color  or  gray  level  coded.  An 
alternative  is  to  color  or  gray  level  code  the  matrix  terrain  data  directly. 
While  analysis  of  these  matrix  image  displays  is  superior  to  trying  to  perform 
analysis  by  visual  inspection  of  the  data  in  printed  numerical  matrix  format, 
they  only  provide  for  a  low  spatial  resolution  analysis  capability.  Shaded 
relief  display  with  variable  illumination  adds  additional  information  for 
analysis  of  all  types  of  matrix  data,  but  is  particularly  meaningful  for  carto¬ 
graphic  data  because  of  the  relationship  to  the  physical  world.  Higher  spatial 
resolution  analysis  of  the  shaded  relief  display  may  be  gained  by  applying 
photogrammetric  models  to  generate  pseudo-stereo-pairs  of  images  (figure  3.)  in 
which  spike  points  are  apparent  under  stereoscopic  analysis.  These  techniques, 
used  singly  or  in  combination,  allow  for  matrix  data  analysis  far  superior  to 
techniques  of  a  decade  ago,  but  they  are  not  enough. 


In  order  to  perform  high  resolution  anomaly  analysis  of  matrix  data  for  the 
purpose  of  either  quality  control  or  information  gathering,  advanced  techniques 
are  required.  These  techniques  include  convolution  filtering,  texture 
discrimination,  and  specialized  color  representation.  n  about  the 
algorithms  used  in  these  techniques  are  described  in  Append! ct w  rough  B. 

* 

Convolution  filters  have  been  used  very  effectively  to  enhance  matrix  data  to 
show  processing  anomalies  as  well  as  where  data  has  been  merged  from  different 
production  equipment,  different  stereo  models,  different  production  methods, 
variable  requirement  specifications,  and  even  from  different  analysts.  These 
types  of  filters  are  used  extensively  by  the  image  processing  community  to  detect 
edge  differences,  and  then  to  re-apply  the  differences  to  sharpen  the  original 
image.  They  also  have  been  shown  to  be  a  powerful  tool  for  the  analysis  of  matrix 
data.  Consider  figures  4a.  and  4b.  Depicted  are  a  shaded  relief  display  and  the 
corresponding  digital  convolution  display  of  a  one  degree  square  of  digital 


FIGURE  SHADED  RELIEF 
OF  TERRAIN  DATA 


FIGURE  4fc.  CONVOLUTION  MASK 
RESULTS  OF  TERRAIN  DATA 
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terrain  elevation  data  produced  from  both  digitized  source  maps  and  stereo-pairs 
of  imagery.  The  shaded  relief  display ,  as  well  as  the  corresponding  gray  level 
coded  elevation  display,  does  not  clearly  suggest  that  the  data  was  produced  by 
two  methods.  In  contrast,  the  convolution  display  clearly  delineates  the 
difference.  In  the  top  third  of  the  display,  the  contouring  effect  due  to 
interpolation  of  digitized  source  maps  is  clearly  visible,  and  is  starkly 
contrasted  with  the  rest  of  the  display,  wherein  stereoplotter  noise  plus  real 
high  frequency  data  combine  to  produce  a  variety  of  edge  differences. 

Texture  analysis  has  proven  to  be  very  powerful  for  discriminating  small  area 
data  that  has  been  merged  into  larger  matrices.  For  example,  digital  elevation 
data  produced  at  a  one  arc  second  interval  was  sampled  and  merged  with  three  arc 
second  data.  This  data  merge  is  not  apparent  in  the  gray  level  coded  elevation 
display  (figure  5a.).  Because  of  inherent  higher  frequency  information  in  the 
one  arc  second  data,  there  is  a  noticeable  difference  between  the  textures  of  the 
levels  of  data  in  the  shaded  relief  display,  and  an  even  more  apparent  difference 
in  the  convolution  display  (see  figures  5b.,  5c.).  A  similar  difference  is 
noticed  when  high  frequency  digital  culture  data  is  merged  with  lower  resolution 
data. 

For  the  purpose  of  determining  compatability  between  data  types,  such  as 
between  digital  terrain  and  culture  data,  simple  color  coding  and  overlay  in  Red- 
Green-Blue  (RGB)  space  may  not  be  sufficient.  A  more  powerful  technique  employs 
coding  each  data  type  along  an  Intensity-Hue-Saturation  (IHS)  axis  and  then 
converting  from  IHS  space  to  RGB  space  prior  to  display.  Since  the  visual 
perception  process  can  distinguish  variation  between  IHS,  the  data  types  can  be 
overlayed  without  a  merging  of  colors,  and  therefore,  without  an  information 
loss.  Consider  figures  6a. -6e.  Figures  6a.  and  6b.  are  gray  level  displays  of 
terrain  shaded  relief  and  digital  culture  information,  respectively.  These  two 
matrices  were  assigned  as  Intensity  and  Hue,  respectively,  with  Saturation  held 
constant.  Figures  6c. -6e.  are  the  Red,  Green,  and  Blue  components  of  the  final 
IHS  color-coded  image. 

Finally,  and  probably  unique  to  map  type  data,  is  the  technique  of  computer 
generating  landform  scenes  as  seen  by  various  visual  and  electro-optical 
sensors.  This  allows  for  a  final  quality  control  analysis  of  information 
content,  and  also  has  been  very  valuable  in  the  definition  of  data  base  require¬ 
ment  specification.  Figures  7a. -7c.  depict  computer  generated  visual 
simulations  and  actual  photography  over  the  data  base  shown  in  the  IHS  figures. 

Impact 

The  impact  of  these  digital  image  technology  concepts  to  matrix  data  display 
and  analysis  lun  been  enormous.  Not  only  has  there  been  a  greatly  increased 
capability  for  the  degree  and  sophistication  of  quality  control,  but  there  is  an 
associated  cost  savings  in  both  the  quality  control  review  process,  and  in  the 
resultant  expense  of  using  matrix  data  containing  anomalies.  These  techniques 
have  application  to  all  types  of  matrix  data  production  as  well  as  for  the 
analysis  of  all  types  of  real  matrix  data. 


i 


FIGURE  6a.  INTENSITY  COMPONENT: 
SHADED  RELIEF  DATA 


FIGURE  6b.  HUE  COMPONENT 
DIGITAL  CULTURE  DATA 


FIGURE  6c .  RED  COMPONENT 


FIGURE  6d.  GREEN  COMPONENT  FIGURE  6e .  BLUE  COMPONENT 


FIGURES  6a-6e* INTENSITY  AND  HUE  INPUT  DATA  TRANSFORMED  TO  RED,  GREEN,  AND 
BLUE  OUTPUT  DATA. 


n&URe  7b.  Visual  Simulation  from  Shipboard  Level  • 


7c.  Actual  Photograph  from  Shipboard  Level 


Figures  Port  Angela*,  Washington  and  Mt.  Olympus:  Comparison 
Between  Visual  Simulations 
and  Actual  Photography 
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Appendix  A.  Elevation  Image  Algorithm 


Digital  terrain  elevation  data  (DTED)  consists  of  terrain  elevation 
values,  in  integer  meters,  stored  on  tape  in  matrix  format.  The  DIED  is 
organized  into  one -degree  by  one -degree  square  geographic  areas  with  a 
maximum  of  ten  one -degree  squares  on  any  one  tape.  For  a  given  area  of 
interest,  not  necessarily  a  one-degree  square,  consider  the  DTED  values 
to  form  an  p  by  q  matrix,  T.  The  first  row  of  such  a  T  matrix,  as  well 
as  all  other  p  by  q  matrices  described  here  and  in  other  sections, 
corresponds  geographically  to  the  southernmost  values  represented  by  the 
matrix.  Consequently,  in  the  schematic  for  a  T  matrix,  as  presented  in 
Figure  A.l,  the  row  convention  is  correct  even  though  it  appears  to  be 
an  inversion  of  the  customary  method  of  displaying  matrix  rows. 


t  t  .  t 

P>1  P>2  p , q 


T 


t  t  . . .  t 

2,1  2,2  2 , q 

t  t  .  t 

1,11,2  1  ,q 


i 

I 

I 


Figure  A.l  Schematic  for  a  T  Matrix 


Let  t  denote  the  DTED  value  in  the  ith  row  and  jth.  column  of  T. 
itj 

Similarly,  let  e  denote  the  corresponding  pixel  m  an  p  by  q  matrix, 

i»  j 

E,  which  is  to  be  recorded  on  film  as  a  digital  image  display  of  T  as  a 
gray  level  coded  elevation  image*  A  gray  level  value  for  each  pixel  of 
E  may  be  determined  by  the  simple  scaling  algorithm  below* 

1.  Compute  a  gray  level  value  for  e  by 


•(G  -  G  )  (A.l) 

max  min 


where : 

t  -  DTED  value  in  the  ith  row  and  jt_h  column  of  T 

i  j  j 

t  =  minimum  DTED  value  in  T 
min 

t  =  maximum  DTED  value  in  T 
max 

G  =*  maximum  gray  level  value  allowed  in  E,  usually  255 
max 

G  =  minimum  gray  level  value  allowed  in  E,  usually  0 
min 

i  53  1,2, ...  ,p 


j  —  l,2,...,q 


2.  Constrain  the  e  from  equation  (A, 1)  to  be  in  the  interval 

i,j 

(G  G  ).  That  is,  set 
min,  max 


(A. 


where : 

G  =  minimum  gray  level  value  allowed  in  E,  usually  0 
min 

G  =  maximum  gray  level  vaLue  allowed  in  E,  usually  255 
max 

i  =  1,2,  .  .  .  ,p 
j  “  1  >  2  ,  *  •  •  }  q 

Note  the  constraints  of  equation  (A. 2)  are  imposed  along  with  a  G  value 

min 

of  zero  to  change  any  e  having  a  computed  value,  on  a  UNIVAC  1100  series 

i,j 

computer,  of  minus  zero  (36-bit  word  of  ones)  to  a  positive  zero  (36-bit 
word  of  zeros).  Further,  this  change  in  e  zero  representation  must  be 

i,  j 

performed  before  the  right -most  8-bit  byte  of  each  e  is  extracted  and 

j 

recorded  on  tape  as  an  output  pixel,  because  extracting  the  right -most  byte 
from  an  e  value  of  minus  zero  would  incorrectly  yield  an  output  pixel 

gray  level  value  of  255  (8-bit  byte  of  ones).  A  similar  rationale  is  the 
motivation  for  the  constraints  imposed  by  equations  (B.3)  and  (C.6)  of 


Appendices  B  and  C,  respectively. 


Appendix  B.  Relief  Image  Algorithm 


Let  r  denote  the  gray  level  value  in  the  ith  row  and  jtji  column 

of  an  p  by  q  matrix,  R,  of  shaded  relief  pixels  over  the  geographic  area 
represented  by  T  (see  Appendix  A).  The  algorithm  below  may  be  used  to 
compute  any  r  provided  T  is  available  and  assuming  the  area  depicted 

j 

by  R  is  illuminated  from  the  Northwest* 


1*  Compute  the  elevation  difference,  ^EL,  between  t  and  it*s 

i,j 

neighbor  t  ,  by 

1HU+1 


ClEL  =  t  -  t  (B.l) 

if  j  i-l,J+l 


where: 

t  =  DIED  value  in  the  ith  row  and  jth_  column  of  T 
i>  J 

t  =  DTED  value  in  the  (i~~l)th  row  and  (j+l)th  column  of  T 

i-l.J+l 

i  =  p  ~1 

53  2,3,#*#,  q-1 


j 


t 


maximum  DTED  value  in  T 


max 

t  =  minimum  DTED  value  in  T 
min 


&g 

min 

N 

g 

f 

b 


i 

j 


minimum  gray  level  step,  computed  as  the  integer  value  of 

(G  -  G  +  1)/(N  -  1) 

max  min  g 

integer  number  of  unique  gray  levels  allowed  in  R, 
usually  33 

round  off  factor,  computed  as  (&g  -  X),  but  set  to 

min  r 

zero  if  the  computed  value  is  negative 

bias  term  used  to  assign  a  background  gray  level  value 

to  any  r  in  a  region  of  relatively  constant  height, 

i  J 

computed  as  the  integer  value  of  (G  -  G  )/2 

max  min 

2.3. ,  p_i 

2.3.. ..,  q-L 


3.  Constrain  the  r 


from  equation  (B.2)  to  be  in  the  interval 


id 

(G  ,G  ).  That  is,  set 
min  max 


r 

G 

r  t  G 

min 

id  min 

< 

►  if  < 

G 

r  >  G 

max 

I 

i ,  j  max 

- 

where : 

G  =  maximum  gray  level  value  allowed  in  R 
max 

G  58  minimum  gray  level  value  allowed  in  R 
min 

i  *  2,3,.*.,  p-1 
j  32  2 , 3,  •  •  ♦  |  q*l 

After  the  r  have  been  determined  by  equation 

id 


(B.3) 

usually  255 
usually  0 

(B  •  3  )  ,  R  may  be 


completed  by  setting  all  elements  in  rows  1  and  p  and  columns  1  and  q 
to  zero  (clear). 


Appendix  C.  Convolution  Image  Algorithm 

Again,  as  discussed  in  Appendix  A,  consider  the  DTED  values  to  form 
an  p  by  q  matrix,  T,  where  t  denotes  the  DTED  value  in  the  ith  row 

i.J 

and  jth  column  of  T*  Similarly,  let  c  denote  the  corresponding  pixel 

in  an  p  by  q  matrix,  C,  which  is  to  be  recorded  on  film  as  a  digital 
image*  Let  H,  an  1  by  l  matrix,  denote  a  filter  to  be  applied  to  T  via 
a  discrete  convolution  process*  Matrix  C  is  produced  by  convolving  T 
with  H  and  subsequently  applying  a  noise -cleaning  and  gray  level  mapping 
process*  Individual  elements  of  C  may  be  determined  by  the  algorithm 
below. 

1*  Form  an  p  by  q  matrix,  T*,  of  elevation  differences  relative  to  the 

minimum  DTED  value,  t  ,  in  T  via 

min 

t’  =  t  -  t  (C.l) 

i,j  l,j  min 

where: 

t  *  DTED  value  in  the  ith  row  and  jth  column  of  T 

t  “  minimum  DTED  value  in  T 

min 


2.  Convolve  T'  with  H  to  produce  an  p  by  q  matrix,  V,  by 


1  L 

v  -!!-<*'  l+1  i+i),(h  } 

i,j  m-1  n=l  m+i-^ySn+j-——  m,n 


(C.2) 


where: 


t' 


2  2 


h 

m,n 


I 


*  elements  of  elevation  difference  matrix,  T1, 
from  equation  (C.l) 

-  elements  of  filtering  matrix,  e.g.,  one  of  those 
listed  in  Appendix  D 

=  an  odd  integer  indicating  the  row  and  column 
dimensions  of  H,  usually  3 


i 


j 


=  i±i 
2  ’ 

=  1±L 
2  1 


+  2, 


+  2, 


P 


1-1 

2 


n  -  ^ 
q  2 


After  the  v  have  been  determined  by  equation  (C.2),  V  may  be 

i,j 

1-1 

completed  by  setting  all  elements  in  the  outer  rows,  i.e.,  1=1,  ...  ,  — 

i  *  p  -  ,  p  and  outer  columns,  i.e.,  j  ■  1,  ...  , 

1-1 

j  *  q  -  i  •••  ,  q  to  zero.  Note  the  total  number,  NQ,  of  elements  in 


the  perimeter  region  of  V  set  to  zero  by  this  process  ace  depicted  as 
the  hachured  portion  of  Figure  C.l  and  may  be  computed  via 


»0  -  a  -  1>*  C(P  +  q)  -  a  -  1)1 


(C.3) 


3.  Apply  a  noise -cleaning  process  to  V  to  produce  an  p  by  q  matrix,  W, 


5.  Constrain  the  c  from  equation  (C.5)  to  be  in  the  interval  (G  ,G  ), 
i,j  min  max 

That  is,  set 


c  =  < 

r 

G 

min 

►  if  < 

f 

C  £  6 

G 

c  >  p 

max 

^  J 

i, j  max 

V  - 

(C.6) 


where: 

G  15  minimum  gray  level  value  allowed  in  C,  usually  0 
min 

G  13  maximum  ^ray  level  value  allowed  in  C,  usually  255 
max 

€  «  threshold  value  chosen  to  force  output  gray  level  values  of 

G  for  those  c  in  a  region  of  DTED  values  of  relatively 
min  i,j 

constant  height,  usually  21 

i  “  1#2,  • • • 9  p 

j  *  •««>  <1 


Appendix  D.  Convolution  Masks 


Twenty-seven  casks  were  selected  for  testing  and  built  into  a 
computer  program,  DMATER,  used  to  generate  test  images  via  these  filters. 
The  first  three  are  Laplacian  edge  detection  masks  while  the  remaining 
twenty -four  masks  allow  for  computing  1st.  through  6th  differences  in 
the  N-S,  E-W,  NW-SE,  and  SW-NE  directions.  Specific  matrix  elements  for 
each  mask  along  with  a  unique  reference  number,  in  parenthesis,  to  its 
left  are  provided  below. 


The  three  Laplacian  masks  are: 

0-10 

(1)  -1  4  -1 

0-10 

-1  -1  -1 

(2)  -1  8  -1 

-1  -1  -1 

1  -2  1 

(3)  -2  4  -2 

1-2  1 


The  masks  for  1st.  through  6th  differences  in  the  N-S  direction  are: 

0-10 

(4)  0  10 

0  0  0 

0  10 

(5)  0-20 
0  10 

0  0-100 

0  0  3  0  0 

(6)  0  0  -3  0  0 

0  0  10  0 

0  0  0  0  0 

0  0  10  0 

0  0  -4  0  0 

(7)  0  0  6  0  0 

0  0  -4  0  0 

0  0  10  0 


-10  0 

(16)  010 

0  0  0 

10  0 

(17)  0-20 

0  0  1 

-1  0  0  0  0 

0  3  0  0  0 

(18)  0  0  -3  0  0 

0  0  0  1  0 

0  0  0  0  0 

1  0  0  0  0 

0  -4  0  0  0 

(19)  0  0  6  0  0 

0  0  0  -4  0 

0  0  0  0  1 

-l  0  0  0  0  0  0 

0  5  0  0  0  0  0 

0  0  -10  0  0  0  0 

(20)  0  0  0  10  0  0  0 

0  0  0  0  -5  0  0 

0  0  0  0  0  1  0 

0  0  0  0  0  0  0 

1  0  0  0  0  0  0 

0  -6  0  0  0  0  0 

0  0  15  0  0  0  0 

(21)  0  0  0  -20  0  0  0 

0  0  0  0  15  0  0 

0  0  0  0  0  -6  0 

0  0  0  0  0  0  1 

The  masks  for  1st  through  6th  differences  In  the  Stf-ME  direction  are: 
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£  COLOR  TRANSFORMATION  FOR  CRT  DISPLAYS 

Color  displays  such  as  CRT  devices  use  the  additive  color  principle 
to  generate  their  spectra  of  colors.  Red,  green  and  blue  shades  are  the 
usual  primaries  in  color  CRT  displays.  When  three  types  of  data  are  to 
be  displayed,  the  most  common  (and  easiest)  process  is  to  allow  one  data 
set  to  control  the  red,  one  the  green  and  one  the  blue.  This  process  is 
valid  for  reconstructing  a  color  scene  from  its  red,  green  and  blue 
components,  but  generally  it  is  not  valid  if  the  objective  is  to  maximize 
the  amount  of  information  perceived  from  each  of  the  three  data  sets  by 
the  human  observer.  This  is  due  to  the  nature  of  the  eye  and  the  processing 
of  information  by  the  visual  system  [l]. 

What  is  needed  to  maximize  perception  is  a  set  of  parameters  which 
specify  a  given  color  and  which  are  nearly  independent  of  each  other. 

Without  going  into  the  details  it  can  be  argued  that  intensity,  hue  and 
saturation  are  parameters  which  satisfy  these  conditions. 

Intensity  is  the  apparent  brightness  of  a  display.  Hue  is  the 
predominant  "color" .  Saturation  measures  how  much  a  hue  is  "diluted"  by 
white.  On  the  C. I. E/  chroroaticity  diagram,  which  maps  color  perception, 
saturation  and  hue  can  be  approximately  represented  by  a  polar  coordinate 
system.  Saturation  is  represented  by  the  radial  distance  of  a  color  from 
white  which  lies  at  the  origin  of  the  coordinate  system.  Hue  is  represented 
by  an  angular  measure  around  the  origin  from  an  arbitrary  reference  hue. 
Intensity  would  be  represented  by  a  vector  perpendicular  to  the  hue  - 
saturation  plane.  With  these  three  parameters  any  color  can  be  uniquely 
defined.  This  method  of  defining  color  is  called  the  IHS  (intensity,  hue, 
saturation)  representation.  It  should  be  noted  that  this  representation 
defines  a  cylindrical  coordinate  system  which  will  be  used  later. 


e.  /  HUE-SATURATION  PLANE 
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The  choice  of  three  primaries,  such  as  those  in  a  CRT  display, 
defines  a  triangle  in  the  hue-saturation  plane.  Any  hue  or  saturation 
can  be  achieved  within  the  triangle  by  varying  the  relative  amounts  of 
each  primary.  Unfortunately,  colors  outside  of  the  triangle  cannot  be 
achieved  unless  new  primaries  are  chosen.  Intensity  on  the  display  can 
be  varied  by  changing  the  sum  of  the  amounts  of  the  three  primaries. 

(This  defines  intensity  and  is  based  on  the  work  done  by  Faugeras.  [l]) 

Since  the  maximum  intensity  of  each  primary  is  limited  in  a  CRT  display, 
the  maximum  intensity  of  any  mixture  is  also  limited. 

Both  the  RGB  and  IHS  color  representations  can  uniquely  and  completely 
specify  any  color,  therefore,  there  must  exist  a  transformation  that  will 
take  the  parameters  of  a  color  in  RGB  and  give  the  parameters  in  IHS  and 
vice  versa.  This  transformation  is  now  derived  with  the  restrictions  of 
CRT  displays  in  mind  and  with  some  simplifying  assumptions. 

Negative  values  are  not  possible  for  the  primary  values,  so  colors 
are  restricted  to  be  in  the  positive  octant.  Since  the  intensity  is 
represented  as  being  proportional  to  the  sum  of  the  three  primary  components, 
an  intensity  coordinate  can  be  defined  as  a  vector  starting  at  the  origin 
and  in  a  direction  equidistant  from  each  RGB  axis.  This  vector  is  perpendicular 
to  any  plane  of  constant  intensity  defined  by 

R  +  B  +  G  =  constant  # 

With  the  intensity  vector  defined  we  can  relate  the  cylindrical 
coordinate  system  of  IHS  space  to  the  rectangular  coordinate  system  of 
RGB  space.  The  transformation  is  from  cylindrical  to  rectangular  coordinates 
with  a  rotation. 
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C.L.  GEOMETRY  RELATING  IHS  TO  RGB 


•  f-  f  ^  RELATION  OF  PRIMARIES  TO  CHROMATICITY 


For  derivation  purposes  the  three  RGB  primaries  are  assumed  to  form  an 
equilateral  triangle  in  the  hue  -  saturation  plane.  "White"  contains  an 
equal  amount  of  each  of  the  primaries.  In  RGB  space  each  of  the  primaries 
is  assumed  to  lie  along  the  axis  of  a  three  dimensional  rectangular 
coordinate  system. 


COORDINATE  SYSTEM  DEFINED  BY  PRIMARIES 


For  the  rotation  to  line  up  the  R  axis  with  the  I  axis. 
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For  the  transformation  from  cylindrical  to  rectangular  coordinates. 

Rf  *  I  -co<I<a> 


Gf  *  Sf  cos  Hf  0<S,<® 

B1  »  S’  sin  H1 


Combining  the  two  transformations 

R  =  7?  (I“s'(cos  H'  +  sin  H')> 

G  =  JL  (I+S'(i±^  cos  H'  +  sin  H')) 

B  =  ^  (I+S’  cos  R’  +  sin  H’)) 

For  completeness,  the  inverse  transformation  is 
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So  far,  the  derivation  has  strictly  been  a  mathematical  coordinate 
transformation.  The  next  step  is  to  account  for  the  physical  constraints 
of  the  color  display  on  the  RGB  values,  more  specifically, 

0<R,B,G<1 

In  this  case 

CKI^1//! 

and 

o<s'<r5i 

A 

The  saturation  must  be  proportional  to  the  intensity  for  the  relative 
amounts  of  red,  green  and  blue  to  remain  constant.  To  visualize  the 
reachable  colors  with  this  transformation  I  is  set  to  its  -maximum  and  the 
hue  -  saturation  plane  is  examined.  ****** 


TRANSFORMATION  LIMITATION  ON  ACHIEVABLE  COLORS 


The  triangle  defines  the  total  spectrum  of  achievable  colors  with  the 
chosen  primaries.  The  inner  circle  defines  those  colors  achievable  with 
the  above  transformation. 


All  the  colors  inside  the  triangle  can  be  reached  if  slight  adjustments 
are  made  to  the  hue  and  saturation.  First  we  change  the  reference  hue  so 
that  it  corresponds  to  pure  red 


then  we  alter  the  saturation  for  different  values  of  hue 


gi  m  /Is _ 

2  sin(— -H' ) 

6 

/3S 

m 

2  sin  (11-H') 

2 

_  _^5s _ 

2  sinCi^l-ft' ) 

6 

These  values  of  H*  and  S’  can  be  used  in  the  transformation  equations 
given  previously. 

This  derivation  seems  to  yield  good  results  on  actual  displays  but 
several  areas  of  study  could  still  be  pursued: 

1)  A  more  precise  approach  to  using  Faugeras*  eye  model. 

2)  In  general  three  primaries  do  not  form  an  equilateral  triangle 
on  the  C.I.E.  chromaticity  diagram. 

3)  Hue  and  saturation  do  not  appear  as  polar  coordinates  on  the 
C.I.E.  diagram. 
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